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Abstract. We examine the merger stage of IRAS 03158+4227, one of the most luminous ULIRGs from the 
IRAS 2 Jy sample. Deep optical images are presented along with high-resolution NIR images and optical low- 
resolution spectra. IRAS 03158+4227 is confirmed as a component of an equal-luminosity binary galaxy with a 
nuclei separation of 47 h^^ kpc. A long lopsided tail emanating from the companion, which harbours an active 
nucleus, is the most prominent peculiar feature of the binary. The results of numerical simulations admit the 
interpretation of this structure as a product of the tidal interaction between the two components. If the infrared- 
activity of IRAS 03158+4227 is also dynamically triggered by this process, this would implicate that the ULIRG 
phenomenon is not restricted to the final stage of a binary merger. Alternatively, the system may be a multiple 
merger where the partner(s) has/have already coalesced. 
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1. Introduction 

The characteristics of galaxies that emit a substantial 
amount of their bolometric luminosity in the far infared 
(FIR) has been a m atter of debate since th eir dis covery 
(e.g., Rieke & Low |1972| ; Joseph & Wright |1985| ; Soifer 
et al. 1987). Particular interest is focussed on the class of 
ultra- luminous infrared galaxies (ULIRGs), i.e. galaxies 
with quasar-like infrared luminosities of Lir > IO^^Lq 
(e.g., Sanders & Mirabel f996l Kennicutt f996l Genzel 



the r ange from 50% t o abo ut 100% (e.g., Aurie re et al. 
1996|; Cl ements et al. |l996| ; Murphy et al. |1996| ; Due et 
al. 19971; Borne et al. 2000 ). Murphy et al. have analysed 



et al. |1998| Rigopoulou et al. |199S| ; Scoville et al. pOOOj ) 



Nuclear starbursts and/or AGN activity, dynamically trig- 
gered by gravitational interactions, are thought to be the 
energy sources in the cores of ULIRGs. Absorption and 
re-emission in the FIR is expected to be a consequence of 
thermalization of radiation by large masses of dust grains. 

It has been demonstrated by a large number of stud- 
ies that a high fraction of luminous and ultraluminous 
infrared galaxies show morphological peculiarities, such 
as tidal debris and double nuclei, indicative of gravi- 
taional perturbations. For instance, all ULIRGs of the 
original IRAS bright galaxy survey (BGS; Sanders et al. 



1988a) show indications of strong gravitational interaction 



or merging. Published values for the interacting/merging 
fraction among the ULIRGs from various studies cover 
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near-infrared (NIR) and visual images for 46 luminous in- 
frared galaxies with Ljr > 8.5 10^"'^ Lq from the IRAS 2 Jy 
sample. After combining their sample with the BGS, they 
find that 95% of the galaxies in the combined sample show 
evidence for current or past interactions. There are only 
three ULIRGs in the Murphy et al. sample that do not, to 
the limits of the images given there, show indications for 
interactions. Among them is IRAS 03158+4227, one of the 
most luminous ULIRGs (log Lir/Lq = 12.55, /eo/Zioo — 
0.95). This seems astonishing since the fraction of per- 
turbed systems among the FIR-bright galaxies is known 
to increase with _Lir (e.g., Sanders & Mirabel 1996, and 
references therein). 

On the numerical side, simulations clearly show that 
the major morphological features observed in many pe- 
culiar galaxies are explained as due to tidal forces during 
galaxy encounters (e.g., Toomre & Toomre 1972|; Ba rnes & 
Hernq uist |1992| ; Bekki & Noguchi |l99|; Spoke |1997| ; Mihos 
et al. 1998). Self-consistent models of tidally disturbed 
galaxies indicate large gas concentrations in the centres 
due to s trong and sud den g aseous inflow (Negroponte 
White |l98|; Noguchi |l 99 ij; Barnes & Hernquist |l99l ' 



1996t Mihos & Hernquist |l996D . The models predict that 



merger-driven gas-dynamics and associated star forma- 
tion may result in spectacular starbursts (e.g. Noguchi & 
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Ishibashi 198(:; Mihos & Hernquist 1996), although there 
is considerable uncertainty about the treatment of star 
formation and of the feedback from young stars in the 
simulations. Galaxy merger and the ULIRG phenomenon 
were tied together in the models by Mihos & Hernquist 
(|199§). 

ULIRGs are widely claimed to represent the final 



stages of merging galaxies. Murphy et al.(1996) give pro- 
jected linear separations of less than a few kpc for a large 
fraction of their sample. In the case of IRAS 03158+4227, 
Murphy et al. noticed a nearby, resolved component which 
they described as "not apparently interacting with the 
primary galaxy". IRAS 03158+4227 was therefore not in- 
cluded by these authors in their sample of double-nuclei 
systems and was subsequently considered a single system 
where the nuclei separation must be smaller than the avail- 
able resolution of 0'.'8. The projected linear distance to the 
nucleus of the galaxy mentioned by Murphy et al. amounts 
to 18" corresponding to a projected linear separation of 
47 kpc (throughout this paper we adopt an Einstein-de 
Sitter comology with Hq = 75kms"^ Mpc^^). There are 
only a few 2 Jy ULIRGs with nuclear distances larger than 
10 kpc and only one other system (IRAS 14394+5332) 
with a separation of about 50 kpc. 

In this paper, we present deep optical imaging, high- 
resolution (adaptive optics) NIR imaging, and optical 
spectroscopy of the remarkable ULIRG IRAS03158+4227. 
The observational data (imaging and spectroscopy) are 
described in Section 2. The results are presented in Sect. 3 
and are discussed in the context of simulations in Sect. 4. 
Finally, conclusions are given in Sect. 5. 

2. Observations and data reduction 

Most of the observations presented here were taken at 
the German- Spanish Astronomical Centre on Calar Alto, 
Spain. A summary of relevant data for all observations is 
given m Table|ll 

Table 1. Log of observations. 



filter/ 


^cxp 


instrument 


epoch 


grism 


[s] 






J 


360 


ALFA/Omega-Cass 


Sep 98 


H 


180 


ALFA/Omega-Cass 


Sep 98 


K' 


90 


ALFA/Omega-Cass 


Sep 98 




700 


CAFOS 


Jul 98 


B 


4 560 


CAFOS 


Jul 98,00 


R 


3 700 


CAFOS 


Jul 99,00 


I 


820 


CAFOS 


Jul 00 


B400 


1200 (Gl) 


CAFOS 


Jul 99 


B400 


3 000 (G2) 


CAFOS 


Jul 99 


R 


1200 


TLS Schmidt 


Jan 99 


I 


960 


TLS Schmidt 


Dec 99 



Deep direct optical imaging was performed with the fo- 
cal reducer camera CAFOS at the 2.2 m telescope in the B, 



R, and I band and without filter. CAFOS was equipped 
with a SITe CCD with a scale of 0'/5/pixel. The condi- 
tions were always photometric with a seeing of typically 
about 1". Unfortunately, IRAS 03158+4227 is located only 
about 2' away from the bright foreground star HD 20489 
(V=8.6) and 7'.'4 from the USNO-2 star No.275-02191303 
(B=16.2, R=14.9). Therefore, we took several sets of rela- 
tively short exposures (typically between 100 and 300 s) to 
avoid saturation effects of the CCD. The total integration 
time amounts to 2.8 hours. MIDAS standard algorithms 
were applied for the data reduction. The combined image 
was PSF-deconvolved using the Lucy-Richardson method 
(MIDAS procedure deconvolve/fluey). 

High-resolution imaging in the J, H, and K' bands 
was performed using the adaptive optics system ALFA in 
combination with the Omega-Cass camera (Hippler et al. 
1998 ) at the 3.5m Calar Alto telescope. Omega-Cass uti- 
lizes a 1024x1024 HAWAII detector and was operated at 
the pixel scale of 0'.'08. Since there is no sufficiently bright 
star close to IRAS 03158+4227 which could be used for 
wavefront sensing and with the laser guide star being not 
operational at the time of the observations, we corrected 
the static aberrations of the telescope on a nearby star and 
imaged the target with the deformable mirror "frozen". 
Although this procedure does not yield diffraction-limited 
resolution, it improved the image quality considerably, 
leading to sub-arcsecond resolution. Two adjacent fields 
were observed for deriving the sky frames. During the data 
processing, the images were rebinned in order to enhance 
the signal-to-noise ratio which led to a final pixel scale 
(from the astrometric solution) of 0'/155. After correction 
for fiat field and bad pixels, the images wer e filter ed using 
the wavelet algorithm of Pantin & Starck ( 1996 ) to min- 
imize noise amplification in the subsequent Richardson- 
Lucy deconvolution. The USNO-2 star No.275-02191303 
served as PSF reference for the deconvolution. JHK' pho- 
tometry was derived from the non-deconvolved images and 
tied to the JHK' magnitudes of the PSF star according to 
its entry in the 2MASS Second Incre mental Release Point 
Source Catalog (Cutrie et al. 2000| ). The derived fluxes 
(Fig.||) refer to a synthetic aperture of 7'/75 diameter. The 
photometric error amounts to 0.06 mag. In order to assess 
the separation of a possible double nucleus, the FWHM 
of the images of the ULIRG in the three filters were com- 
pared to those of stars in the field. The average stellar 
FWHMs derived from Gaussian fits amount to 0'/60, 0'.'57, 
and 0'.'58 for J, H, and K', respectively. The FWHMs of 
the ULIRG are 0'.'74, 0'/70, and 0'.'71. This leads to beam- 
deconvolved sizes of 0'/44, 0'/40, and 0'.'40 for the angular 
extent of the emitting core region of the ULIRG. 

Low-resolution spectra of both IRAS 03158+4227 and 
its nearest neighbour galaxy were taken with CAFOS 
equipped with the grism B 400 which is suitable for the 
wavelength range A = 3 200... 8 000 A. With a slit width of 
l'/2 the spectral resolution is about 20 A. 

Finally, we observed the field of IRAS 03158+4227 
with the Tautenburg Schmidt telescope at moderate see- 
ing of about 2". The Schmidt camera was equipped with 
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Fig. 1. IRAS error ellipse of IRAS 03158+4227. 
Coaddition of R- and I-band images taken with the 
Tautenburg Schmidt camera. The size of the image is 
200" X 200"; N is up, E is left. SI to S3 designate stars, 
Gl and G2 indicate galaxies. 

a 2kx2k SITe CCD with pixel size of 24^nix24^m 
which yields a field size of 40'x40'. These images were 
used only to evaluate the large-scale environment of 
IRAS 03158+4227. 

3. Results 

Figure |l| shows a moderately deep optical image of the 
field around IRAS 03158+4227. The data for the IRAS 
error ellipse were taken from the IRAS Point Source 
Catalogue, the coordinates for the objects on the opti- 
cal image are from the astrometric solution for the field of 
the P erseus cluster of galaxies (Brunzendorf & Meusinger 
1999 ). The IRAS error ellipse clearly overlaps with the 
optical image of a faint galaxy (Gl) with z = 0.13 
(NASA Extragalactic Database, NED). At the distance 
of the galaxy, the size of the field shown in Fig. ^ corre- 
sponds to about 0.5 MpcxO.5 Mpc. There is no sign for 
a dense galaxy cluster or a rich group of galaxies around 
IRAS 0315+4227. Several small and faint galaxy images 
are seen in the field, but a substantial fraction may be 
dwarf galaxies of the Perseus cluster in the foreground. 

Murphy et al. (1996) presented a K-band image taken 
with the Palomar 200 inch telescope at a seeing of O'.'S 
that does neither show evidence for a neighbour galaxy 
outshined by the brighter foreground star SI on a less 
resolved image, nor for a double nucleus of Gl. As men- 
tioned in Sect. 1, IRAS 03158+4227 was assumed to be a 
late merger with nuclei separation of le ss tha n O'.'S, even 
though it was noted by Murphy et al. (199£) that there 




Fig. 2. Spectral energy distribution (SED) of 
IRAS 03158+4227. Asterisks denote our measure- 
ments while diamonds mark the IRAS fluxes. (The data 
point at 12 /im represents an upper limit.) 



is a galaxy (G2) at a distance of 18". The two galaxies 
Gl and G2 have approximately the same apparent mag- 
nitudes; after calibrating our B-band image using stars 



from the USNO-A2.0 Catalogue (Monet et al. |1998D we 
find B = 18.6 and 18.4 for Gl and G2, respectively. This 
corresponds to M-q k, —21 A and —21.6 (for Ab = 0.8 mag 
from the NED||_ and a fc-correction fee = 0.5 mag from 



Coleman et al. 198C) and to a very strong FIR excess of 
logLpiR/iB = 2.3 for Gl (see e.g., Soifer et al. 1987 for 
the definitions of Lfir and Lb). 

The spectral energy distribution (SED) of 
IRAS 03158+4227 is shown in Fig.|. The SED is 
typical for ULIRGs, with the bulk of the luminosity radi- 
ated at far-infrared wavelengths. In addition to the data 
points shown in Fig.|, we note that IRAS 03158+4227 
is identified with a radio continuum source in the NVSS 
(Condon et al. [19981). The flux density of 12.4mJy at 



1.4 GHz corresponds to a ra dio-to -infrared fiux ratio (for 
definition see Helou et al. 1985) q = 2.64. This is in 



agreement with the well-known correlation between the 
fiux densities in the infrar ed and the radio continuum (see 
Sanders & Mirabel 1996 , and references therein) where 
{q) « 2.35 for most of the luminous infrared galaxies 
in the BGS, but is larger for the galaxies with highest 
infrared luminosities. (For instance, Helou et al. found 
q = 2.60 for Arp220.) 

The optical images (Fig. ^ clearly reveal that G2 has a 
faint, but very extended material arm on the side opposite 
to Gl. This feature is seen in all optical bands, and the 
Lucy- Richardson-deconvolved image shows that it is quite 
narrow as expected for tidal tails. In the optical bands, 
the arm is about 3 mag fainter than the main body of G2. 
The surface brightness of the brightest part of the tail is 



^ The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and 
Space Administration. 
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Fig. 3. Optical images of IRAS 03158+4227: a B-band; b R-band; c I-band; d co-addition of all optical images; e the 
same image as d, but after Lucy-Richardson deconvolution with 10 iteration steps; f the contour plot of image e where 
successive contour lines correspond to a factor 3 in intensity. (The unfiltered image shows the same structures as the 
images a to c and is not displayed here). The size of each image is 30" x 42", N is up, E is left. 



estimated to /ib ~ 25 mag arcsec"^, the projected linear 
extent is about 70kpc. Illustrative examples for similar 
tidal structures with one dominating arm can be found 
among the Arp ( |1966D and Arp & Madore ( |1987D galax- 
ies (e.g., Arp 107, 110, 129, 173, 252, 255, AM 2350-302, 
AM 0552-324, AM 0427-475). For the host galaxy of the 
ULIRG itself, tidal signatures are much weaker, although 
there appears to be some fuzz around Gl. Several brighter 
spots may represent tidal debris, perhaps including huge 
star forming regions; but clearly, none of these spots is 
bright enough to be considered a second nucleus of Gl. 
Unfortunately, it can not be excluded that faint morpho- 



logical features of Gl are hidden behind the bright star 
SI. 

The radial intensity profiles of Gl and G2 in the R- 
band are shown in Fig. ^. Gl may be approximated by two 
components with different scalelengths. Alternatively, it 
can be classified as an "E-like" ULIRG: the radial surface 
density profile is reasonably fit by a deVaucouleurs r^/^- 
law over the range R . . .30kpc. (The innermost 1"5 
were e xclude d from the analysis.) According to Sanders 
et al. ( ^000|) , about one third of the ULIRGs from the 
1 Jy sample are classified as "E-like". Contrary to Gl, the 



Meusinger, Stecklum, Theis, Brunzendorf: IRAS 03158+4227 



5 






1 1 1 1 

J. 






¥ 




M -1 

O 






-2 


III 






10 20 


30 



-1 



1 


1 1 


- i 






T 










1 


S 

1 1 



1.6 



r (kpc) 



-2 



-2 



1.8 2 2.2 
[r (kpc)]'-'^ 



2.4 



10 15 
r (kpc) 



[r (kpc)]'' 



Fig. 4. Intensity profiles of Gl (top) and G2 (bottom) 
derived from the R-band image by means of the MIDAS 
procedure fit/eU3. Iq is the central intensity. The angular 
interval where the image of Gl is affected by the star SI 
was excluded from the profile analysis. 



the image of SI and a further nucleus, but the probability 
for such a configuration is very low. 

The optical spectra of Gl and G2 (Fig.^) show strong 
emission lines; the equivalent widths (EWs) are listed in 
Table H A radial velocity difference between Gl and G2 
of 240kms~^ is estimated, but the uncertainty due to the 
low spectral resolution is quite large (about ±200kms~^). 

IRAS 03158+4227 has sometimes been classified as a 
Seyfert2, whereby strong absorption was invoked as the 
reason for the absence of indications of nuclear activity 



in the hard X rays (Risaliti et al. 200C). In the spectrum 
from Fig.^, we do not see clear-cut evidence for an AGN 
in IRAS 03158+4227; in particular the conventional diag- 
nostic EW ratio [O iii]5007/H/3 is too low for a Seyfert 
nucleus. However, the signal-to-noise ratio is low for these 
lines, and, owing to the low dispersion, we are not able 
to resolve Hq:+[N II]6584. The low intensity of the emis- 
sion lines in the blue part of the spectrum is probably 
best explained as due to strong internal dust absorption 
and resultant reddening, though the effect of the under- 
lying stellar absorption is difficult to estimate. If H/3 is 
significantly affected by stellar absorption of an older, A- 
type burst population, the intrinsic [O iii]5007/II/3 ratio 
would be even less compatible with a Seyfert spectrum. 
On the other hand, the diagnostic line ratios of G2 are 
more consistent with an AGN. It should be noted that the 
K'-band image of Gl is stronger concetrated than that of 
G2 (Fig.||) as seems typical for ULIRGs as compared with 
AGN hosts (Sanders et al. ^OOOj ). 



Table 2. Equivalent widths EW (in units of A) of the 
emission lines measured for the galaxies Gl and G2. 









Gl 


G2 


EW(H/3) 


3.5 


± 


0.4 


2 + 2 


EW([Oiii]5007) 


8.0 


± 


0.4 


59 + 1 


EW([0 i]6300) 


8.0 


± 


0.4 


0.1 + 0.2 


EW(Ha+[Nll]6584) 


118.7 


± 


0.2 


124.0 + 1.0 


EW([Sii]6717,6731) 


3.7 


± 


0.3 


1 + 1 



profile of G2 is better approximated by an exponential law 
which suggests a "disk-like" structure. 

The images in the NIR-bands are shown in Fig.||. 
These high-resolution images were taken to search for 
point sources indicating either a close double (or mul- 
tiple) nucleus with small separations or a nucleus out- 
shined by the bright star SI in the optical image, rather 
than to evaluate extended structures of low-surface bright- 
ness. For such an aim, observations at longer wavelengths 
are needed, since the morphology of the central parts of 
ULIRGs is strongly affected by dust obscuration. As ex- 
pected (cf. Sanders et al. 2000| ), the light distribution in 
K' is very compact for the galaxy Gl which appears as a 
point source. Figure|5|does not reveal any sign for a double 
nucleus in Gl or for a nucleus close to SI down to the res- 
olution (beam-deconvolved size) of 0'.'40. Of course, there 
remains the possibility of an incidental superposition of 



4. Discussion 

The long, curved tail of G2 seen in Fig. ^, along with the 
fact that Gl and G2 have nearly the same redshift, may 
be taken as an indication for the gravitational interac- 
tion of Gl and G2. Such tails of escaping debris from 
the far side of a victim disk are well known indicators 
of the encounter of nearly equal-mass spiral galaxies (e.g. 
Toomre & Toomre |1972t Schombert et al. |1990| ). On the 
other hand, if the ULIRG activity was triggered by this in- 
teraction, the large projected distan ce bet ween these two 
galaxies is surprising. Murphy et al. ( 1996| ) suggested that 
ULIRGs with large (> 10 kpc) nuclear separation may rep- 
resent triple merger with a third, undetected nucleus from 
a previous encounter or, alternatively, that the ULIRG 
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Fig. 5. High-resolution images of IRAS 03158+4228 in the J, H, and K' bands. The scale and the size are the same as 
in Fig.|. 



phenomenon can occur in an early phase of interaction. 
Below, we briefly discuss IRAS 03158+5228 in the light of 
these two scenarios. 

Arguments in favour of the multiple merger scenario 
were derived from the properties of elliptical galaxies 
(e.g., Bar nes |1984| Mamon |l98^; Schweizer |1989| ; Weil & 
Hernquist 1996), from the dynamical diversity of ULIRGs 
(Borne et al. pOOO| Cui et al. |200lD , and from de tailed 
studies of individual galaxies (Taniguchi & Shioya 1998; 
Lipari et al. 2000). Multiple encounter and merger are 
suggested t o occu r natura lly in compact gro ups of galax- 



ies (B arnes |1989| ; Hickson |l997t Borne et al. |2000t Bekki 
2001 ). It seems likely that a fraction o f ULI RGs is trig- 
gere d by such a process. Borne et al. (2000) and Cui et 
al. (2001) considered the appearance of double or mul- 
tiple nuclei as a key-test for the multiple merger origin 
and derived percentages of 20% and 17%, respectively, of 
multi- nuclei ULIRGs. The fraction of ULIRGs triggered 
by multiple merger is certainly larger than the fraction of 
multi-nuclei systems since a multiple nucleus is expected 
to evolve on a short timescale to a double nucleus and 
finally to a single nucleus. Unfortunately, this method is 
faced with serious difficulties which can lead to an over- 
estimation of multi-nuclei systems: the morphology of the 
central regions of ULIRGs has the tendency to be strongly 
affected by dust obscuration effects and by the appearance 
of regions of intense star formation on a scale of kpc or sub- 
kpc. Further, the studies mentioned above did not identify 
real interacting members with spectroscopic observations. 



Following Bekki (2001), it seems fair to say that the frac- 
tion of ULIRGs formed by multiple merging is still highly 
uncertain. 

Dinh-V-Trung et al. ( ^00 1]) have studied the SIX sys- 
tems with nuclear separations larger than 20 kpc among 



the ULIRGs from the complete 1 Jy sample (Kim & 
Sanders 1998| ). The optical and K'-band imaging obser- 
vations and optical spectra suggest the multiple merger 
scenario only for one of those ULIRGs, IRAS 14394+5332. 
It can not be excluded that IRAS 03158+4227 is a mul- 
tiple merger hke IRAS 14394+5332. Indeed, the mor- 
phologies of these two systems show some similarities. 
Though we do not find evidence for a close double nu- 
cleus, IRAS 03158+4227 might be in a more advanced 
stage, where two nuclei of Gl have coalesced and the in- 
ner region is already well relaxed as seems indicated by 
the radial luminosity profile of (Fig.|^). 

How can the nuclear activity of G2 be matched by such 
a scenario? Nuclear activity is not unusual in ULIRG- 
systems, though the active nuclei are mostly located in 
the hosts of the ULIRG itself. However, there is one sys- 
tem (IRAS 17028+5817) among the widely separated pairs 
studied by Dinh-V-Trung et al. where the spectrum of the 
ULIRG's host is of Hll-type whereas the companion has 
a LINER-type spectrum. We can not exclude that G2 is 
also a late merger. However, its disk-like structure admits 
a variant of a multiple merger where the ULIRG activity 
was triggered by a past merger and the AGN in G2 by the 
present interaction between Gl and G2. 

The simulations by Bekki ( ^001 ) have demonstrated 
that multiple merger can trigger repetitive starbursts with 
a star formation rate comparable to ULIRGs. However, 
the discussion by Bekki suggests that very intense star- 
bursts with an amplitude of ^ lO^Mo/yr are not likely 
in such an environment. As was stressed already in the 
Introduction, IRAS 03158+4227 is one of the most lumi- 
nous ULIRGs from the 2Jy sample. According to the 
relation derived by Clements et al. ( |l996|) , the 60 /im 
flux transforms into a huge star formation rate of about 
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Fig. 6. Optical low-dispersion spectrum (not 
flux-calibrated, observer frame) of the ULIRG 
IRAS03158-t-4227 (Gl, top) and of the galaxy G2 
(bottom). 



As an alternative to the multiple merger scenario, it 
seems tempting to speculate that the activities in the cen- 
tres of Gl and G2 were triggered by the same process, 
namely an interaction of Gl and G2. Liu & Kennicutt 
( |1995| , their Fig. 4) discussed the empirical distribution 
of the equivalent widths of the Ha-|-[Nll] line for differ- 
ent merger morphological types. The EWs measured for 
Gl and G2 (Table |) are in better agreement with Liu & 
Kennicutt's morphology type 3 (= systems of two disk 
galaxies) than with type (2 — advanced merger) which 
appear to be single. Moreover, according to its infrared 
colour index /25//6O < 0.2, IRAS 03158-^4228 belongs 
to the group of "cool" ULIRGs which are characterized 
as major merger with prominent extended tidal struc- 
tures and resolved double nuclei rather than by sma ll 
(< 2.5 kpc) nuclei separation systems (Surace et al. 20001 ). 

The simulations by Mihos & Hernquist ( 1996 ) have 
demonstrated that disk/bulge/halo systems with dense 
central bulges experience strongest gaseous infall and star 
formation activity in the final stages of coalescence when 



they are within a few kpc of one another. Their disk/halo 
models without dense bulges, on the other hand, are most 
active in earlier phases of merging when the galaxies are 
separated by tens of kpc. At the beginning of the first star- 
burst phase, the snapshots of the disk/halo merger mod- 
els by Mihos & Hernquist (their Figs. 11 and 12) show 
a remarkable similarity with the few morphological de- 
tails seen in IRAS 03158-1-4227: one galaxy (hereafter: gl) 
is more concentrated, especially the gas and the young 
stars, with knots and short arms whereas the most promi- 
nent feature of its interaction partner (hereafter: g2) is 
an extended curved tail at the opposite side. During the 
next time steps, when the SFR reaches its maximum, the 
bridge between gl and g2 becomes weaker and gl becomes 
more concentrated. 

The long lopsided tail of G2 is the only visible mor- 
phologically peculiar feature of the system. It is therefore 
important for the understanding of the merger stage of 
IRAS 031584-4227 to know whether this structure can be 
due to the tidal interaction with Gl. Since it is not pos- 
sible to follow the evolution of the extended tidal struc- 
tures in the snapshots shown by Mihos & Hernquist we 
performed a small series of restricted N-body simulations 
like those in Toomre & Toomre (1972). The main idea of 
this method is to derive the orbits of both galaxies from 
the corresponding two-body problem, e.g. by solving a 
Kepler problem, if the galaxies are treated as point masses. 
Using these orbits the time-dependent potential at each 
point is given by a superposition of the two galactic po- 
tentials. Stars are treated as test particles, by this reduc- 
ing the classical iV-body problem to N single-body prob- 
lems. ( Detai ls of the applied code are described in Theis & 
Kohle (2001)). The main advantages of this method are a 
fast computation and a high spatial resolution. However, 
the method is not self-consistent, because effects of self- 
gravity (like fragmentation or dynamical friction) are ne- 
glected. Anyway, comparisons between self-consistent and 
restricted N-body calculations demonstrated in several 
cases a good agreement, provided the encounters are not 
too strong and/or the duration of the simulated stage is 
not very long. Therefore, and because there are not many 
constraints from observations, the restricted A^-body sim- 
ulations should be a good starting point. The results from 
the present simulations arc considered however indicative 
rather than conclusive. 

Figure shows the result of a parabolic encounter with 
an orbit inclination of 60°, a minimum distance of 15 kpc 
(reached at t = 0), a line-of-sight-velocity of 200 kms~^, 
and a final projected distance of 50 kpc. The total (dynam- 
ical) mass of each galaxy is IO^^Mq. Since we are mainly 
interested in the tidal shape of G2, we resolved only G2 in 
test particles. Motivated by the observed radial intensity 
profile, we distributed the test particles in an exponential 
disk with a scale length of 5 kpc and a cut-off radius at 
15 kpc. For simplicity we assumed the disk to be coplanar 
with the plane of sky. The qualitative agreement between 
the optical image and the simulations is very good: af- 
ter 7.5 • lO'^yr a lopsided structure has been formed. This 
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Fig. 7. Simulation of the gravitational perturbation of the galaxy G2 by the galaxy Gl. Both galaxies have the same 
dynamical mass of 10^'^ M©. Left: projected distribution of the mass particles from G2. Gl is marked by the cross at 
the bottom. The length is measured in kpc and the time unit is 1.49 Myr. No test particles have been used inside the 
central kpc. Right: blurred intensity plot of the image from the left hand side, overlaid by a contour plot; the lowest 
intensity contour approximately corresponds to the detection threshold in the optical images of IRAS 03158-f 4227. 



gives rise to an arc-like structure north-west of G2 and a 
small spur emanating north-east from the main body of 
G2. In the simulations, these features are the dense parts 
of a tidal tail which seems also to exist in the observations 
(e.g. lower-right panel of Fig. ^). The bridge connecting Gl 
and G2 may be too weak to be clearly detected. Increasing 
the line-of-sight velocity to 300 kms~^ yields very similar 
results. 

The assumption that IRAS 03158-1-4227 is triggered 
by the interaction between Gl and G2 implicates that 
the ULIRG phenomenon is not restricted to late binary 
merger stages. Such an interpretation is supported by 
further indications. Rigopoulou et al. ( 1999| ) reported a 
lack of any correlation between the stage of merger, mea- 
sured by the separation of nuclei, and the infrared lu- 
minosity in an unbiased sample of 62 ULIRGs. Further, 
there is no trend for increased ULIRG activity in systems 
with more centrally concentrated Ha emission (Mihos & 
Bothun 1998), and also the total mass of molecular gas 
in ULIRGs is not related to the linear separation (Gao & 
Solomon |1999|; Ri gopoulou et al. |1999D . Finally, Dinh-V- 
Trung et al. ( pOOl] ) present evidence for IRAS 23327-^2913 
to be hosted by a non-disturbed spiral-like galaxy which 
may be interpreted as an early stage of merging. 



equal-luminosity giant galaxies separated by 47 hj^ kpc. 
The companion galaxy, G2, harbours an active nucleus 
and has an extended, curved tail. This tail may be inter- 
preted as a tidal structure induced by the gravitational 
interaction with the host galaxy of the ULIRG. The high- 
resolution NIR images do not reveal any sign for a close 
double nucleus down to a (resolution-limited) nuclei sep- 
aration of 1 kpc. 

ULIRGs in such widely separated systems may be ex- 
plained by a multiple merger, like IRAS 14394-f 532 (Dinh- 
V-Trung et al. 2001), where the interaction of two compo- 
nents has already reached an advanced stage with a very 
small separation of nuclei. Indeed, the optical luminosity 
profile of Gl is reasonably fit by a de Vaucouleurs law. The 
profile of G2 is more compatible with a disk-like structure. 
Hence it seems likely, in the multi-merger scenario, that 
the ULIRG in Gl was triggered by a past merger and the 
AGN in G2 by the present interaction between Gl and G2. 
If the nuclear activity of G2 is also due to a (late) merger, 
the system IRAS 03158-1-4227 would represent the case of 
the gravitational interaction of two advanced merger. Such 
scenarios would implicate that multiple merger in compact 
groups can trigger very intense starbursts. 



5. Conclusions 

Deep optical images, spectra and high-resolution NIR im- 
ages are presented for IRAS 03158+4227, one of the most 
luminous ULIRGs of the 2 Jy sample. The host galaxy, 
Gl, is identified with a component of a binary of nearly 



As an alternative explanation, both the ULIRG 
IRAS 03158+4227 and the nuclear activity in its neigh- 
bour galaxy might be dynamically triggered by the same 
process, namely the gravitational interaction of these two 
galaxies. In this case, IRAS 03158+4227 would represent 
an early stage of binary tidal interaction, and the inter- 
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pretation of ULIRGs as final merger stages may need to 
be reexamined. 
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